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Abstract 

Recent experiments by the CLEO III detector at CESR indicate that the J/ip spectrum produced 
in T decay is in conflict with Non-Relativistic QCD (NRQCD) calculations. The measured J/ip 
momentum distribution is much softer than predicted by the color-octet mechanisms. The expected 
peak at the kinematic limit is not observed in the data. However it has recent been pointed out 
that NRQCD calculations break down near the kinematic endpoint due to large perturbative and 
non-perturbative corrections. In this paper we combine NRQCD with soft collinear effective theory 
to study the color-octet contribution to the T — > J/ip + X decay in this region of phase space. We 
obtain a spectrum that is significantly softened when including the correct degrees of freedom in 
the endpoint region, giving better agreement with the data than previous predictions. 
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I. INTRODUCTION 



Since the discovery of the J/i/j, heavy quarkonium decay and production have served as 
a laboratory for testing both perturbative and nonperturbative aspects of QCD dynamics. 
The large mass tuq of the heavy quark inside the quarkonium set a high mass scale (~ 2ttt,q) 
at which the QCD approaches the asymptotically free regime and thus one could hope to 
calculate perturbatively. On the other hand, hadronization process happens at the much 
smaller mass scale of order rriQV 2 , where v is the typical velocity of the heavy quarks in 
the quarkonium. For quarkonium, rriQV 2 is numerically of order Aqcd so the involved QCD 
process also contains nonperturbative effects. 

During the past 15 years, the interactions of the non-relativistic heavy quarks inside the 
quarkonium have been understood to some extent using the framework of Non-Relativistic 
QCD (NRQCD) [J 3- NRQCD is an effective theory that reproduces the physics of full 
QCD by adding local interactions that systematically incorporate relativistic corrections 
through any given order in the heavy quark velocity v. This effective theory provides a 
generalized factorization theorem that includes nonperturbative corrections to the color- 
singlet model, including color-octet decay mechanisms. All infrared divergences can be 
factored into nonperturbative matrix elements, so that infrared safe calculations of inclusive 
decay rates are possible jjjj]. NRQCD solves some important phenomenological problems 
in quarkonium physics. For instance, it provides the most convincing explanation to the 
surplus J/tj) and ip' production at the Tevatron in which a gluon fragments into a color- 
octet cc pair in a pointlike color-octet S-wave state which evolves nonperturbatively into the 
charmonium states plus light hadrons. The NRQCD factorization formalism allows these 
fragmentation procedures to be factored into the product of short distance coefficients and 
long distance NRQCD matrix elements among which the leading one is (0^^,^[ 3 Si\) where 
are l° ca l 4-fermon operators in terms of the fields of NRQCD. 

There are, however, some problems that remain to be solved. One challenging problem is 
with the polarization of J/ip at the Tevatron. The same mechanism that produces the J /if) 
described above predicts the J/i[> should become transversely polarized as the transverse 
momentum p± becomes much larger than 2m c |5|. Though the theoretical prediction is 
consistent with the experimental data at intermediate p±, at the largest measured values of 
p± the J/ip is observed to be slightly longitudinally polarized and discrepancies at the 3a 
level are seen in both prompt J/ip and ip' polarization measurements jsj. 

Recently, a new problem arose as a result of measurements of the spectrum of J/ip 
produced in the T(1S) decay by the CLEO III detector at CESR Q. NRQCD calcu- 
lations have been made for the production of J/ip through both color-singlet and color- 
octet configurations (!, |||. Theoretical calculations predict that the color-singlet process 
T(IS') — > J/ipccg + X features a soft momentum spectrum. Meanwhile, the theoretical 
estimates based on color-octet contribution indicates that the momentum spectrum peaks 
near the kinematic endpoint j9[. In contrast to the theoretical predictoins, the experimen- 
tally measured momentum spectrum is significantly softer than predicted by the color-octet 
model and somewhat softer than the color-singlet case 0. 

The NRQCD predictions break down in the endpoint region because the effective field 
theory does not contain the correct degrees of freedom to describe the physics. NRQCD 
contains soft quarks and gluons, but it does not contain quarks and gluons moving collinearly. 
The correct effective theory to use in situations where there is both soft and collinear physics 



is Soft-Collinear Effective Theory (SCET) [10Hl3j. A combination of SCET and NRQCD 
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has been successful in reproducing the shape of the measured J / if) momentum spectrum in 
e + e~ — > J/if>+X 14|. SCET has the power to describe the endpoint regime by including the 
light energetic degrees of freedom. In addition, renormalization group equations of SCET 
can be used to resum large perturbative logarithmatic correctoins. Nonperturbative martix 
will occur naturely in deriving the factorizatoin theorem using SCET. 

In this paper, we use SCET to study the color-octet contribution to the T — > J /if) + X 
decay near the endpoint. We derive the factorization theorm in SCET for this process. We 
find that the spectrum is significantly softened when including perturbative up to leading 
logarithms (LL) and nonperturbative corrections near the endpoint, giving better agreement 
with the data than the previous predictions. 



II. FACTORIZATION AND MATCHING 

In this section, we briefly derive the SCET factorization theorem for T — > J /if) + X near 
the endpoint. A more detailed derivation will be presented in the Appendix. The derivation 
is similar to radiative T decay |16| , which we refer to for details. However, for the process we 
discussed here, it involves the decay of a heavy quarkonium into another heavy quarkonium, 
thus we should combine SCET with two independent NRQCD's for these two onia systmes. 
The factorization for a similar process B — >■ J /if) + X s has been discussed in Ref. [151 ] . 

Near the endpoint regime, a new factorization formula is required since the NRQCD does 
not include all the relevant physical degrees of freedom and thus the factorization theorem 
breaks down. This can easily be seen when we analyze the kinematics at the endpoint. To 
do so, we work in the centre-of-mass (COM) frame, and introduce light cone coordinates. 
By introducing the parameter x = (E^ + p^)/My, we have 

n My My „ 

Pr = ~Y n ^ + + K , 
Ml „ xMi 



u _ — v a , — l r -ft h n 
2xM T U + 2 n 



Px 



M 



T 



(l- r -^Jn»+ (l-aOn" + fc£ . (1) 



Here n = (1, 0, 0, 1) and n = (1, 0, 0, —1), we have defined r = ml/ml, and we also assumed 
that = 2m c and My = 2m^. and are the residual momentum of the QQ pair inside 
the T and J/if> respectively. Near the kinematic endpoint, the variable x — > 1 and thus the 
jet invariant mass approaches zero. In NRQCD, an expansion of k^/rrix is performed and 
hence the jet mode is integrated out, which is only valid when it has a large invariant mass, 
i.e., away from the endpoint. As x — > 1, the jet becomes energy is large, but the invariant 
mass becomes small, with k^/mx of order 1. Hence we must keep k^/mx to all orders. 
As a result, the standard NRQCD factorization breaks down at the endpoint. SCET is the 
appropriate framework for properly including the collinear modes needed in the endpoint in 
order to make reasonable predictions. 

To derive the factorization theorem in SCET, we start from the optical theorem in which 
the decay rate can be written as 

2E ^=Y^]^J2 J d%e-^(T\O\y)\J/iP + X}(J/iP + X\O(0)\T) 1 (2) 
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FIG. 1: QCD production amplitude for T — > J/ip + X. The J /if) is produced in a color-octet and 
becomes a color-singlet by emitting a soft gluon. There is another contribution to this process with 
only one gluon emitted, which is suppressed by an order of a s . 



where the summation includes integration over phase space of X. The SCET operator O is 



O 



(3) 



where the Wilson coefficient Cj is obtained by matching from QCD to SCET at some hard 
scale n = hh and the SCET current function Jiiu) is contrained by the gauge invariance. 
For instance, in our case, to leading order the non-vanshing SCET current will be of the 
form 



abc 



(4) 



Here F^f" / is a hard coefficient containing the color and spin structures which is obtained 



by matching onto the QCD Feynman diagrams shown in fig. [TJ The matching gives 



abc 



2ig 4 1 1 



•V, 1 - r.\l T Ml 



(5) 



where we have chosen the coefficient so that the Wilson coeffcient C(/i, u) is 1 at the hard 
scale fin- The A's boost the J /if} or T from the COM frame to a frame where those quarkonia 
have arbitrary four-momentum. if> and x are the heavy quark and antiquark fields which 
create or annihilate the constituent heavy (anti-) quarks inside the quarkonia. The collinear 
gauge invariant field strength is built out of the collinear gaugle field A% 



—WUK + 9s(A%, q )i_)W n 
9s 



(6) 
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where 



w n= Yl exp \ 9s ~p n ' An >i ) ' ( 7 ) 



perms 



is the collinear Wilson line. The operator V picks out the large momentum label 12| . 

We can decouple the usoft modes from the collinear degrees of freedom by making the 
field redefinition [13J 

K, q -> YK«Y ] , (8) 

where Y is the usoft Wilson line made out of the usoft gauge fields. In such a way, we can 
separate the collinear physics from the usoft and obtain the factorization theorem in SCET. 
Meanwhile using similar arguments presented in Ref. [l5[ allows us to factorize the matrix 
element into the convolution of the shape functions for T in color-singlet configuration 
and J/ip in color-octet one. Following this procedure, the decay rates can be written as a 
convolution of soft shape functions and the jet function with an overall hard coefficient. By 
introducing z = E^/m b , we get the decay rate of the form 

^ = T P[z,r]J2\C(u^)\ 2 J dk + J dl + Uk + )S4l + )S T (M r (l-x)-k + -l + ), (9) 



where P[z, r] = 8ny/z 2 — 4r/(l — r) is a kinematic factor and 



i^B^^^iw^. (10) 



We have used spin symmetry 17 



AfAf(...a\..^...) = 

^A^' (...a k ...a k ...), (11) 

to simplify the matrix elements, and applied the identity 5^AfA^' = (v s v s ' — g 65 '), where v s 
is the four- velocity of the T or J/ip. 

The shape function for J /if) is defined as 

r dv - il+ (0| \xhiYYT k Y^ Vcdr) 4^ 4*iYYT k YW X A |0> 

where we have made the field redefinition in Eq. OH]) for the two collinear gluons in the final 
state by introducing two different usoft Wilson lines Y and Y. For T we have the shape 
function 

Sril) -J^ e "ww^FT () 

respectively. Both shape functions are normalized so that J dZ + S , ^ x(^ + ) = 1- 
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FIG. 2: Feynman diagram for the leading-order jet function. The spring with a line through it 
represents a collinear gluon. 



The jet function is given by 

(0| [B a a X Bf{ y )B^B h ^)] |0> 

= J (g a a>gw5 aa '5 bb ' +g a p>g ( 5 a >5 ab '5 ba ')5 l 



^ 2 \y^)5(y + )e- L 2 k+ y-Uk + ). (14) 



To the leading order, the jet function can be calculated by evaluating the diagram shown in 
Fig. [2] , which gives 



in 



-0(Af 



(15) 



where Q is the total four momentum carried by the jets. 
The J/ip shape function can formally be written as 



Sf(l 



X. 



\ ai YYT k YW i)J{in ■ d - l + )ala^^laiYYT k Y^ Xc] |0) 



4m c (0«[3^]) 



(16) 



and to lowest order in v 2 , Sr(l + ) S(l + ). By integrating over k + and / + in Eq. dHJ), we find 
the tree level decay rates become 



1 T-l 

— = T P[z,r}Q(l-x) 



(17) 



with P[z,r] = P[z,r]/87r. This can easily be seen to reproduce the tree level calculation of 
NRQCD 0. 



III. RUNNING 

Effective field theories provide a powerful tool to sum logarithms by using the renormal- 
ization group equations (RGEs). Logarithms of the ratio of different scales arise naturally in 
perturbation theory, which can cause a breakdown of the perturbative expansion when those 
scales are well separated. By matching onto an effective theory, the large scale is removed 
to be replaced by a running scale \x and the effective operators are run from a high scale to 
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FIG. 3: One-loop order diagrams needed to calculate the counterterm to the color-octet opera- 
tor. The double line presents the J/ip fields in color-octet configuration while the spring lines are 
the soft gluons. 

the low scale using the RGEs, which sum all large logarithms of the ratio of scales into an 
overall factor. 

In our case, there are logarithms of the form log(l— x) that will appear in the perturbation 
series. Near the endpoint, x — > 1, these become large, and must be resummed. In this 
section, we will apply the RGEs of SCET to sum these large logarithms. 

In the previous section, we have matched QCD onto the SCET operator by intergrating 
out the hard scale fin, replacing it with a running scale fi. We now run the operator from 
this hard scale to the collinear scale. To do so, we calculate the counterterm for the operator 
to determine the anomalous dimension, and then use this in the RGEs. 

The one-loop corrections to the SCET operator in Eq. (Tj0) is given by the graphs in Fig. [31 
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Evaluating these diagrams gives the divergent term 



A- 



loop 



E 



u s C A 
1 

+ - 

e 



11, 

-T + - 2 + lOE 



n-Q 2 jr 



e e 

Co>(n-0 + cu) , n-Q + w u(n-Q-u) , n-Q — u 

— log — — — log — 

n-Q(n-Q — u) 2n-Q n-Q(n-Q + u) 2n-Q 



x Aq. 
(18) 



The calculation lets us estimate the hard scale be fin = n-Q/ \fr which will minimize 
the logarithm. The divergent piece must be canceled by Z 8 Z^/Zo — 1, where Zq is the 
couterterm for the operator in SCET, Z 3 is the gluon wave function counterterm 



1 + (%c, 

47r e 



4n F 



3 3 

and Z% is the counterterm of color-octet J ftp operator 

, , a sCA 



(19) 



Arte 



(20) 



This leads to 



'0 



E 



4tc 



+ 



4 + - (log 



n-Q 2 jr 



1 [U 4n F T F 



uj(n-Q + u) , n-Q + u u(n-Q-u) , n-Q — u 

■ — log — — — log — 

n-Q{n-Q — u) 2n-Q n-Q(n-Q + co) 2n-Q 



(21) 



From Eq. ( I2ip . we can extract the anomalous dimensioin of the operator through the stan- 
dard method. Using the anomalous dimension in the RGE for the color-octet Wilson coef- 
ficient and running from the hard scale down to the collinear scale gives 



|C(£,/i c 



2C A 





b 






a s (n 2 c ) 


n-Q 2 jr 




_a s (n-Q 2 /r)_ 




_a s (n-Q 2 /r)_ 



Hi] 



where 



14 4n F T F 

y ~ ~vc~ A 



-(2e-i) 



log(l - o 



log£ 



(22) 



(23) 



with £ = {n-Q + u)/(2n-Q), b = 11Ca/3 — 2n F /3, and the collinear scale fi 2 c ~ m 2 x . 

At the collinear scale, the jet mode can be regarded as large and can be integrated out. 
The decay rate can be further run down to the soft scale /i s . To do this, we first note that 
the decay rate can be modified to 



d/ + d/ + '6(M T (l -x)- l + ) U S {1 + - 1+', // c , fi s )S4l + ', fi s ) 



(24) 



s 






FIG. 4: One loop corrections to the J/ifi soft shape function defined in Eq. ([12 



since in Ref. 18 , it was shown that 



(r v>6 e(m .d + M r (i-x)- 1+) 4a iX - b \r) = 

0(M T (1 -x)- l + )(T\xla^ b 4a iX - b \T) . (25) 



And here we introduced an evolution kernel U s as in Ref. [19|. The soft shape function has 
evolution through U s which will sum the large logarithms between /i s and \x c . 

The evolution kernel can be calculated explicitly [19], once we figure out the anomalous 
dimension for the soft shape function. To one-loop order, calculating the diagrams in fig. HI 
we get 



Ztt 



1 1 n 
-— - - log 



rM\ 



Therefore we hi, is 



u s {i + -i + ') 



fi s T(-u) 



- 5 (k + ) + - — 
el y 1 e M T 



M T 9(A;+) \ 
k+ I 



(26) 



(1+ - 1+'] 



l+OJ 



(27) 



with uj is defined as 



2C A r s(Mc) ada 4C Al a s (u c 

■log- ' 



where /3[a s 



7 as(Ms) /9[a] &o 
■(11CU/3 - 2n F /3)a 2 /(27r). Note that w < 0. We defined a function K. 



(28) 



k ~ n -Ufa) P[<x] 



[1 + logr) 



2CV a s (/x 
— — ( 1 + log r ) lo 



(29) 



which is related to K in Eq. ( I27j) by 

2C A P s ^ e) ada r da' 



K = 



17 6ga a (// c ) l« s (/i 
The soft scale is /i 2 . ~ rM|(l — x) 2 . 



L s ( Ms ) £H ^ S (M S ) /V] 

txC a ( a s (n c ) _ a s (fi c ) j ^(/i c ) 

a s (/i s ) « s (/i s ) 



(30) 
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Gathering all the pieces we have, we find the resumed decay rate 



dr 

dz 



— = T P[z,r] 



' ii s r(-u) 





2C A 

h) 




8ttC A 
as (KQ2/r)b% I* 1 

/ de 

Jo 






n-Q 2 /r 




_a s (n-Q 2 /r)_ 


_a s (n-Q 2 /r)_ 





dZ + d/ + '6(M T (l -x) -l + ) 



(/+ -i+'y+* 



S4l + ',fi s ). (31) 



IV. PHENOMENOLOGY 

The decay rate from the previous section, Eq. (13T|) . summed up the leading logarithmic 
corrections which are important near the kinematic endpoint. Away from that region, the 
logarithms that we have summed are not important and contributions that we neglected in 
the endpoint become important. We therefore would like to interpolate between the leading 
order color-octet contribution in NRQCD away from the endpoint and the resummed result 
near the endpoint. We thus define the interpolated differential rate as 

dr .„ N /dr\ /dr\ 

ay V dy / nrqcd Vdy/ SCET 

Here, in order to compare with the data, we have used the scaled momentum defined as 
V — Pip/p™ ax - We see that as y — > 1 the first term vanishes, leaving only the SCET 
contribution in the endpoint region. 

To proceed, we need the soft shape fucntion of J/ip that appears in Eq. ( 13~TT) . We will 



apply a modified version of a model used in the decay of B mesons [20 



f(l + ) = {f^(v- l) ab ~ 1 e-<»-Mv - 1) , (33) 

with r) = l + /A. Here A = — M is of order Aqcd; and a and b are adjustable parameters 
of order 1. In our case, we choose a = 1 and b = 2. A was determined so that the first and 
the second moments of the shape function 

POO 

mi = / df + /(f + ) = A(6+l), 

J A 

m 2 = £° dl+(i+f /([+) = A 2 Q + (b + l) 2 ^) , (34) 

take the value 890MeV and (985MeV) 2 respectively. 

We show the results of resumming in fig. [5j The short dashed line is the NRQCD decay 
rate only and the dotted line is the NRQCD decay rate convoluted with the shape function. 
The thin line includes only the perturbative resummed interpolated decay rates without 
convoluted with the soft shape function. The solid thick line presenting our final result is 
the interpolated decay rate convoluted with the shape function in Eq. (1331) . As can be seen, 
the shape function and the perturbative resummation both result in a softer spectrum. The 
combination of the two is softer still. 

In fig. El we compare our results with the experimental data from CLEO Q. We use 
the values m c = 1.5 GeV, m& = 4.9 GeV, and Aq CD = 0.2 GeV so that a s {2m h ) = 0.1793. 
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lAW/dy 
1.0 r 




FIG. 5: Comparison between the NRQCD results and SCET predictions normalized to the NRQCD 
decay rate at the end-point. Here y = Pip/p™ ax is the scaled momentum. The short dashed line 
is the NRQCD decay rate only and the dotted line is the NRQCD decay rate convoluted with the 
shape function. The solid thin line includes only the perturbative resummed interpolated decay 
rates without convoluted with the soft shape function. The solid thick line presents the interpolated 
decay rates convoluted with the shape function. 



The solid line represents the color-octet interpolated decay rate convoluted with the shape 
function. The NRQCD matrix element was chosen to be (X | f 1 ^] | T) = 2.3GeV 3 . For 
comparison, we have used in the plot the same value for the overall strong coupling evaluated 
at the scale 2m c as in Ref. j9[. The shaded band is obtained by varying the NRQCD 
color-octet matrix element from 0.003 GeV 3 to 0.014 GeV 3 . Since the numerical value 
of the matrix element (C^[ 3 Si]) is fixed by experimental data, it has large uncertainties 
coming both from experiments and theoretical higher order corrections. For comparison, we 
also show the color-singlet contribution as the dashed line 0, ls[ . The complete spectrum 
involves a combination of the color-octet contribution we calculated here, and the color- 
singlet component @, Q shown in the figure. 

The differential rate predicted by the color-octet model is peaked near the end-point 
region. When convoluted with the shape function, the momentum distribution is shifted to 
the left but still close to the kinematic limit. Once we resum the large leading logarithms 
under the framework of SCET, interpolate the result with the NRQCD prediction and then 
convolute with the soft shape function, we find that spectrum is significantly softened near 
the endpoint and the peak is pushed further to the left, in better agreement with the data. 
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FIG. 6: Comparison of the color-octet contribution to the differential rate to the data from 
CLEO [7j. The solid thick line presents the interpolated decay rates convoluted with the shape 



function with a choice of (C?[ 3 5i]) = 6.6 x 10 3 GeV 3 [21[. The shaded band is obtaiend by varying 



M 3 Si]) from 0.003GeV 3 to 0.014GeV 3 [2|. Here, we also show the color-singlet contribution in 



long dashed line 0, S|- The complete spectrum will involve a combination of both the color-octet 
and color-singlet contributions. 



We note here that if we use a high scale for the overall coupling constant, the color-octet 
contribution will be much smaller. In order to make a consistent comparison of theory 
to data, one needs to treat the endpoint of the color-singlet contribution in SCET and 
NRQCD, which we have not done here and will leave for future work. In the perturbative 
expansion, the color-singlet process is suppressed relative to the color-octet one by an factor 
of a s . However, there is a large enhancement due to kinematic factors in the diagram, which 
can be as large as 360 that provides a huge enhancement to compensate the perturbative 
suppression [8[]. 



V. CONCLUSION 

In this paper we studied the color-octet contribution to T — > J/ip + X near the kinematic 
endpoint. In this regime, the usual NRQCD factorization formalism breaks down due to large 
perturbative and nonperturbative corrections. We combined the usual NRQCD effective 
theory with SCET, which contains the correct physical degrees of freedom, to derive a 
factorizatioin theorem for the differential decay rate, dT/dy, valid in the endpoint region. 
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This also allows us to resum large logarithms which appear in the endpoint by running the 
rate from the hard scale to the collinear scale and then to the soft scale using the RGE 
of the effective theory. At the soft scale, we are left with NRQCD shape functions. Using 
models for the color-octet shape function, and interpolating away from the endpoint to the 
leading order NRQCD prediction, we are able to make predictions about the color-octet 
contribution to the decay rate over the entire kinematic region. 

Though a quantitative comparison to data can be made only when including both the 
color-singlet and color-octet terms, some qulitative conclusions can still be drawn here. 
We note that once we sum the large leading logarithms using the framework of SCET 
and convolute with the soft shape function, the spectrum is significantly softened near the 
endpoint and the peak is broadened and shifted to the left. This effect greatly improves the 
agreement between the data constraints and the theoretical predictions. 

We note here that the hard scale we have chosen in our final results is 2m c the same as in 
[9| . This scale is much smaller than what we expected from our one loop order corrections. 
If we used the hard scale \xu decided by our caculations, the color-octet contribution will 
be further suppressed. As a result, we expected that the dominant contribution for this 
process will be from color-singlet rate. Once including higher order effects as well as the 
feed-down of i/j(2S) and Xcj to J /if) [7J, the color-singlet decay rate should also be softened. 
We expected that combining color-singlet and octet contributions will give a good fit to the 
data. 
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Appendix A: Deriving The Factorization Theorem for the Endpoint 

In this Appendix we show the factorization theorem for T — > J/ip + X near the endpoint. 
We begin with the differential decay rate 

2E *^ = T^M^E / d'ye-^ (T\0\y)m + X) (J/if> + X\O{0)\T) , (Al) 

where X includes both the usoft X u sector and the collinear X c sector. The operator O is 
the SCET operator defined in Section [Til of the form 

Q = e -i(M r v+-pf).y V a^u jabpT ^3^] Q^S^ , (A2) 

where 

rrab nnl r>b± ( \ o\ 

a T [l 3 Si] = ^(Ai-^xs, (A4) 
0?[8 8 fiy = ^(A 2 -a) M T c Xc -. (A5) 
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Inserting the operator into Eq. ( lAlj) . the 0^{y) picks up an additional phase and the 
differential rate becomes 

- 16^M r V abc a ' h ' d 

x (ti j$ a Tt [i 3 ^i] <y*[&si\{ v )\jii> + x) 

W + X\ J^ l [l 3 S 1 } O^[8 3 ^](0)|T) 

pOs/3/ii/t -pa'/3'fi'u' Aabc,a'b'c' / \c?\ 

= L abc 1 a'b'c' ^ot^p,v,a' >' V ' V AD / 

In the exponent of Eq. ( 1A6I) . we have used g M — M T f M + Pn fJ '/2 »s M T /2(1 — As 
mentioned in Section HH we can decouple the usoft modes from the collinear degrees of 
freedom using the field redefinition [l3| 

An, q YA^ q Y , (A7) 

which modifies C^[8 3 5i] to 

of ->■ yyo^yW = ojf . (A8) 

Using this field redefinition, we can write 

A abc,a'b'c' 1 f H 4 -iM r /2(l-x)n-y 

^."'^c'"' ~ 16tt 3 M t 7 27 

x (T| a Tt [l 3 Si] 6^[8 s 5 a ](y) «* a* O^l 3 ^] a^[8 3 ^](0) |T> 

x(0|^(y)^;(0)|0), (A9) 

where we have introduced an interpolating field, a^, for the J/^ and used the completeness 
of states in the usoft and collinear sectors 

(AlO) 



J2\X C )(X C \ = 1. (All) 



Xr 



The the usoft Wilson lines only come with the color-octet J/ip operator, 0°f[8 3 S\). The T 
is a very compact bound state, due to the large 6-quark mass. In a multipole expansion, 
long wavelength gluons interacts with the T color charge distribution through its color 
dipole moment since the state itself is color neutral. In the theoretical limit of very heavy 
bottom quark, this coupling to the dipole vanishes 15j. The order of the corrections can be 
estimated by means of the "vacuum-saturation approximation" [l| . A complete set of light- 
hadronic states can be inserted between the T operator and the J/ip operator. 

Notice that the T operator is in color-singlet configuration, therefore the sum over states is 
saturated by the QCD vacuum |0) with corrections of order u 4 [l[. 
Thus we are able to write 

A abc,a'b'c' _ 1 f J4 p ~iM Y /2(l-x)n-y 

~ lQ n 3 Mr J a V e 

x (T\O^[l 3 S l ](y)Ol[l 3 S 1 ](0)\T) 
x(0\d^[8 3 S 1 ](y)a^d^[8 3 S 1 ]\0) 
x<0|J#(y)J$(0)|0>. (A12) 
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Now we substitute the expressions for the T and J/ip shape functions as well as the 
jet function from Section [Til into Eq. (1A12I) . and use pidp^/(2E^,) = mly/z — Ardz/2. The 
result for differential decay rate becomes 

^ = T P[z,r] ^|CV,/i)| 2 I dA; + I d/ + J w (k + ) S^(l + )S r (M r (l - x) - k + - /+)(A13) 
Thus we obtain the desired factorization theorem in the endpoint region. 
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